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Distortions in the primordial cosmic microwave background (CMB) along the line-of-sight can be 
modeled and described using 11 fields. These distortion fields correspond to various cosmological 
signals such as weak gravitational lensing of the CMB by large-scale structure, screening from 
patchy reionization, rotation of the plane of polarization due to magnetic fields or parity violating 
physics. Various instrumental systematics such as gain fluctuations, pixel rotation, differential gain, 
pointing, differential ellipticity are also described by the same distortion model. All these distortions 
produce B-mode that contaminate the primordial tensor B-modes signal. In this paper we show 
that apart from generating B-modes, each distortion uniquely couples different modes (li 7^ I2) 
of the CMB anisotropies, generating {EB) correlations which for the primordial CMB are zero. 
We describe and implement unbiased minimum variance quadratic estimators which using the off 
diagonal correlations in the CMB can extract the map of distortions. We perform Monte-Carlo 
simulations to characterize the estimators and illustrate the level of distortions that can be detected 
with current and future experiments. The estimators can be used to look for cosmological signals, 
or to check for any residual systematics in the data. As a specific example of primordial tensor 
B-mode diagnostics we compare the level of minimum detectable distortions using our method with 
maximum allowed distortion level for the B-modes detection. We show that for any experiment, 
the distortions will be detected at high significance using correlations before they would show up 
as spurious B-modes in the power spectrum. In another words if B-modes were to be detected, the 
presence of any non-primordial contribution could be checked and removed using the off-diagonal 
correlations in the CMB. 

I. INTRODUCTION 

Inflation is perhaps the most promising paradigm of the early universe. Inflation resolves various problems with 
the standard Big-Bang scenario (see [ll, Q for reviews) and provides seed scalar perturbations for structure formation. 
Besides the scalar perturbations, inflation also predicts tensor perturbations or gravitational waves. Detection of these 
tensor perturbations is considered as a "smoking gun" for inflationary scenarios. Primordial scalar perturbations create 
only E-modes of the CMB^, while primordial tensor perturbations generate both parity even E-modes and parity odd 
B-modes polarization [1, @, 0] ■ The detection of primordial tensor B-modes ^ in the CMB would confirm the existence 
of tensor perturbations in the early universe. Various observational efforts are underway to detect such B-mode signal 
of the CMB gil. The magnitude of tensor generated B-mode signal is directly related to the Hubble parameter H 
during inflation, and thus a detection would establish the energy scale at which inflation happened. The amplitude 
of primordial B-modes can be characterized as a ratio of tensor-to-scalar perturbation amplitudes, r. The limit from 
WMAP 5-year data is r < 0.22 

Apart from foreground subtraction challenges, and the challenge of reaching the instrumental sensitivity to detect 
primordial B-modes, there are also many sources of non-primordial B-modes which contaminate the inflationar y sig nal. 
Some example physical processes that generate B-modes are, weak lensing of the CMB by large-scale structure 



11^. 



patchy re- ionization [isl.fl^. rotation of the CMB polarization [l^,|T^. Instrumental systematics also generate spurious 
B-modes [13 • What if B-modes are detected? How would we convince ourselves that the B-modes are primordial? 
In this paper we show that any line-of-sight "distortion" in the primordial CMB that generates non-primordial B- 
modes also leaves a distinct signature in the otherwise zero {EB) and {TB) correlations d, 0, 0, [3 ■ We show that 
the amplitude of {EB) (and TB) correlations is proportional to the distorting field while the B-mode power spectrum 
is proportional to the power spectrum of the distorting field. We construct an unbiased minimum variance quadratic 
estimator for the distortion field T) using the quadratic combination of observed E, B and T CMB maps. We find 
that if an experiment sees non-primordial B-modcs generated by distortion T), our estimator would detect it at higher 
significance. 



^ To first order in perturbations, primordial scalar perturbations do not generate B-modes of CMB. However at second (and higher) order 
in perturbations, scalar perturbations do produce B-modes The B-modes generated from higher order perturbations are expected 

to be smaller than the tensor B-mode levels that the upcoming and future experiments (like CMBPol) are sensitive too. 

^ Throughout the paper the term primordial B-modes refers to the B-modcs generated during inflation via tensor perturbations. 
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The organization of the paper is as follows: In Sec. |TT] we introduce a model for describing the distortions in 
the CMB, and discuss the effect of distortions on the CMB fields and their correlations. In Sec. IIIII we construct 
and characterize a minimum variance estimator which can reconstruct the distortion fields. In Sec. IIVI we show the 
correspondence between our 11 distortion fields with the known physical and instrumental effects on the CMB. In 
Sec. |V] we discuss numerical results and discuss the use of our formalism, with an emphasis on primordial B-mode 
diagnostics. Finally in Sec. IVII we summarize our findings. 

II. DISTORTIONS IN THE PRIMORDIAL CMB MAPS 

We parametrize the distortions in the CMB polarization map following Ref . [l3| . Changes in the Stokes parameter 
of the CMB due to distortions along the line-of-sight can be written as 

6[Q ± lU] (n) - [a ± i2lu] (n) [Q ± iU] (n) + [/i ± ih] (n) [Q T lU] (n) + [71 ± il2] (n)f (n) 

+crp(n) • V[Q ± tU]{n; a) + <j[di ± zd2](n)[9i ± id2]f{n; a) + a\{n)[di ± id2]'^f{n; a) + . . . . (1) 

The first line captures the distortions in a single perfectly known direction n. The distortions in second line capture 
mixing of the polarization fields in a local region of length scale a around n. We Taylor expand the CMB fields Q, U, 
and T around the point n and consider the leading order terms. Here Q, U, and T stands for primordial (un-distortcd) 
CMB fields. Since {Q ztiU){n) is spin ±2 field, a(n) is a scalar field that describes modulation in the amplitude of the 
fields in a given direction n; ix'(fi) is also a scalar field that describes the rotation of the plane of polarization, (/i ±1/2) 
are spin ±4 fields that describe the coupling between two spin states (spin-flip), and (71 ± 172 )(n) arc spin ±2 fields 
that describe leakage from the temperature to polarization (monopole leakage hereon). Distortions in the second line 
of Eqn. dll), {pi ±^2), (i^i ± ^2), and q are measured in the units of the length scale a. The field {pi ± ip2){n) is a spin 
±1 field and describes the change in the photon direction; we will refer to it as a deflection field. Finally {di ± 'i2)(n) 
and q{n) describe leakage from temperature to polarization. (c?i ± c?2)(n) is spin ±1 field and we will refer to it as 
dipole leakage; g(n) is a scalar field that we will call quadrupole leakage. 

These 11 distortion fields correspond to various cosmological distortions such as weak gravitational Icnsing of the 
CMB, screening effects from patchy reionization, rotation of the plane of polarization due to magnetic fields or 
parity violating physics and various instrumental systcmatics such as gain fluctuations, pixel rotation, differential 
gain, pointing, differential ellipticity are also described by our distortion model. In Sec. IIVI wc will discuss the 
correspondence between the distortion flelds and instrumental and cosmological signals in more detail. 

We will work in the fiat-sky limit where scalar fields such as the CMB temperature T and a complex field (iSi±iiS2)(n) 
of spin ±s can be expanded in the Fourier basis as 

T(l) = J dnT(n)e-*' ", 

[5i±i52](l) = (±l)^ydn[5i(n)±i52(n)]eT^*^'e-''-", (2) 

where (p\ ~ cos^^(n • 1). The complex field {Q ± iU){h) is a spin ±2 field, whose Fourier harmonics are referred as 
(E ± iB){l). We want to calculate the change in the CMB field due to the distortion field 'D{h). 

We Taylor expand the observed CMB fields assuming that the distortions are small, keeping terms only to leading 
order in the distortion field V. The changes in the E and B fields due to the modulation a, rotation oj, spin-flip 
(/i, /2) and deflection (pi,P2) take the following form 

S(L) = / ^V{h)E{h)W^{h,h), 

E{L) = E{L)+ I -flLv{l,)E{h)W^ih,h), (3) 

where I2 = L — li and , W-^ are given in Table HI We have assumed zero primordial B-modes for our flducial 
model. Similarly the changes due to temperature leakage, (71, 72), 9, and (^1,^2) can be described by 

BiL) = I |^I?(ll)f(l2X(ll,l2), 

E{L) = ^(L)+ J ^V{h)f{h)W^{h,h). (4) 
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TABLE L Filters, f'^^'i'^uh) 



V 


.fEB{h,h) 




/?b(1i,12) 


W^#(li,l2) 


W^#(ll,l2) 


a 


C',^^sm2(v5ii -¥^1, 


) 


C,^^sin2(v5ii -<^i,) 


sin[2(<^i2 - i^l)] 


C0S[2(VP12 - <y9L)] 


LO 


2(7,-^^ cos 2((y3ii - ifii 




2C','^^cos2(v9i, -<^i,) 


2cos[2(v3i2 - v^l)] 


-2sm[2(v3i2 - V^l)] 


71 


Cf^'^ sm2(</3L -ipi2 


) 


C'iT^ sin2(</PL - "y^ij 


Sm[2(93li - <y3L)] , 


cos[2(</?ii - (/jl)] 


72 


Cl^'^ cos2(<y3L - V'la 


) 


C";^^ cos2(v3L - ^12) 


cos[2(v3ii - , 


-sin[2(v3ii - (^l)] 


/i 


C',^^sin2(2(^L - Vii - 




Cf^^ sin2(2(pL - '^ii - ¥^12) 


sm[2(2(^ii - ipi^ - (/pl)] 


cos[2(2i^ii - V3i2 - v'l)] 


h 


C'/;^^cos2(2^L -^1, - 




C','^^cos2(2(pL -¥^11 -V'ls) 


cos2(2iy3ij - V3i2 - V'l) 


-sin2(2iy3ij - v'la - Y'l) 


di 


Cfi'^ihcr) cos((^L + V'li - 


2^1 J 


Cf^ {ha) cos{ipi^ + ipi^ -2ifi^) 


-(l2(j) COs[(/3ij + ipi^ - 2v3i] 


-(bcr) sin[(/3ij + (/3i2 - 2i^l] 


d2 


-Cf'^^iha) sin((/3L + ■y^ii ^ 




~ClJ'{lia) sin(v3L + Vii ~ 2(^12) 


(Z20-) sm[(^ij + V3i2 - 2i^l] 


(/2O") COs[(/?lj + 1,312 ^ 2(^l] 


q 


-C',^^(Zia)^sm2(v3ii - 


^12) 


-C,7(/ia)2sin2(^,, -^,,) 




— (^20-)^ cos[2((^i2 - </?l)] 


pi 


-(7,-^^(7(11 X L)sm2(v3ii 


~ 'PI2) 


-^^^^(li X L)sm2(<y3ii -lySla) 


a{h X ii) ■ z sm[2(i^i2 - (^l) 


a(l2 ■ii)sm[2(<^i2 -Vl)] 


P2 


-Q^^a(li -L) sin2(v5i, 


- VI2) 


-(5(^^(7(11 ■ L)sm2(v3ii - (^ij 


cr(l2 • li) sin[2(((9i2 - 


cr(l2 X li) • zsm[2(v5i2 - 



It is clear from Eqs. ([3]) and (H]) that the distortion field of wavevector L mixes the polarization modes of wavevectors 
li and I2 = L — li. 

The power spectra of the different modes of the primordial CMB field are defined as 

= {27rrSih+h)Cf^^' , (5) 

where X, X' = T,E,B and (E{li)B{h)) = 0, (f (li)B(l2)) = . The off diagonal terms of the power spectrum {h ^ I2) 
are zero because of statistical isotropy and the (£'(li)_B(l2)) and (T'(li)_B(l2)) correlations are zero because is 
parity even, is parity odd, and the physical processes responsible for the generation of the CMB anisotropics are 
parity conserving. The observed CMB power spectra in the presence of distortion and lensing effects can be written 
as 

(X(li)X'(l2)) = (2^)2<5(li+l2)[C^^'+AA,f^"], (6) 

with 

where a is the FWHM of the beam and and Ap are the temperature and polarization detector noise, respectively. 
In our numerical calculations we assumed fully polarized detectors, for which Ap = \/2At. Using Eqn. ([3]) and (|4]) 
we can calculate (X(li)X'(l2)) correlations. Taking the ensemble average over the CMB fields for a fixed distortion 
field 2?(n), one gets 

(X(li)X'(l2)>CMB - l2)I?(ll + I2) , (8) 

for XX' = TT, TE, TB, EB and 

r / 7^Q^''C',^''^#(1M")<(1',1") P = a,a;,(/i,/2),(pi,P2), 

(S(ll)i?(l2))cMB,P = (271)2^(11 + l2)<^ ^ (9) 



/ 7TrT2 1") 25 = temperature leakage, 
J (27r) 



(2^) 

where 1" = li — 1' and we have assumed zero primordial B-modes for our fiducial model. The filters fxx'O-iih) are 
given in Table[Tl Note that (...)cmb means an average over CMB realizations. The leading order effect of distortions 
on the {E{li)B{l2)) and (r(li)i3(l2)) correlations is to introduce nonzero off-diagonal terms, while for {B{li)B{li)) 
the leading order effect is to generate diagonal terms. 
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III. MINIMUM VARIANCE ESTIMATORS FOR DISTORTION FIELDS 

Our goal is to use Eqn. ([8]) to construct a suitable estimator for the Fourier components 2?(L) of the distortion 
field in terms of the observed fields T{l),E{l), and B{1). Following wc can define a minimum variance quadratic 
estimator Vxx' (L) for I?(L) by weighting quadratic combinations of different CMB fields X by window functions 



Vxx'iJ^) = A^^^,{L) J ^X{h)X'ih)F^x'ih,h), 



where L = I2 + li , and the normalization 

A'^x'iL)^ 



J (27r)^ fxx'0-i^h)Fxx'ihA2) 



(10) 



(11) 



is chosen to make the estimator unbiased, i.e. (X'(L))cmb — ^(L). The weights F^^, can be optimized by minimizing 
the variance subjected to the normalization constraint 



d 



dF-^^.ihM) 



Vxx' - V 



0. 



(12) 



CMB, distortion 



The angle bracket stands for average over both the CMB and the distortion fields. For XX' = TT, EE and BB, the 
minimization yields 



F^xih.h) = 



For XX' = TB and EB, for which Cf^' = 0, we get 



F^x'ih,h) 



f^xih,h) 



f^X'ih,h) 

cxxrX'X' ' 



(13) 



(14) 



and for XX' = TE we find 



Fte{^iM) 



/t£;(1i' I2) 



r'TTr'EE I r'TEr'TE > 



where Cf^-^ and are the observed power spectra including the effects of both the signal and the noise. 

The variance of the estimator is 



(I?XX'(L)I?yy'(L')) CMB, distortion — (2.)^^(L + L'){ 



XX', YY 



>{L)+N^^'^,'yy,{L)+N^ 



XX'. YY 



where 



^^XX',YY' 



(2^ 



(15) 



(i) + ...}, 



(16) 



Fl^,(h,h) F?y,{\,M)Cf:^ CfJ +F?y,{hA,)Cf^^CfJ ,(17) 



with L = li + I2. For the minimum variance estimator, the Gaussian noise N'^^, xx'i-^) = -^xx'iF) ~ A^-x,{L) 
provides the dominant contribution to the variance. The higher order noise terms X^^j^Tyy, (L) come from the 
connected part of the trispectrum and are proportional to (C'^'^(i))". To a good approximation these non-Gaussian 
noise terms are much smaller than the Gaussian noise N'^-^,{L). In Fig. ([T]) we show the Gaussian noise N'^g{L) for 
all the 11 distortion fields, we discuss these plots in detail in Sec.|Vl Although for the experiments like CMBPol, EB 
estimator is the most sensitive, for experiments such as PLANCK which has lower sensitivity to polarization, one can 
combine all the estimators (TT, TE', EE, EB, TB) to construct a minimum variance estimator. 
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Signal-to-Noise: To calculate the signal-to-noise we need to assume a form for the power spectrum of the distortion 
field C^^{L). Motivated by the discussion in [i3j[1Hj we assume that the distortion fields, as defined in Eqn. ([1]), 
arc statistically isotropic and Gaussian. Thus their statistical properties are fully described by their power spectra, 



(r>(l)P(l')> = (27r)2j(l + . 
For illustration we assume power spectra of the form 



•DP 



exp(-Z(; + l)a^/81n2). 



(18) 



(19) 



i.e. white noise above a certain coherence scale ap. The parameter Ax> characterizes the rms of the distortion field 
v. The signal-to-noise for a given distortion field can be calculated as 



Signal 
Noise 



a 



vv 



N, 



vv 



1/2 



(20) 



where Cf^ is the distortion power spectrum for a distortion T) given by Eqn. (HH), N'^'^{i) is the Gaussian noise for 
the estimator of the distortion V (given by Eqn. ((T7])). In Fig. [2] we show total signal-to- noise for all the 11 distortions 
as a function of imax- In Scc.|V]we discuss the plot in detail. 



IV. CORRESPONDENCE BETWEEN THE DISTORTION FIELDS AND COSMOLOGICAL AND 

INSTRUMENTAL EFFECTS 

A. Cosmological Signals 

The distortions defined by Eqn. ((T|) are general distortions along the line-of-sight of the CMB. Here we consider 
known physical processes which generate such distortions in the CMB. 

Cosmological rotation: The plane of linear polarization of CMB field can be rotated due to interactions which 
introduce a different dispersion relation for the left and right circularly polarized modes, during propagation from the 
surface of last-scattering to the Earth. Such interactions can come from three main sources: (a) interactions with 
dust foregrounds, (b) Faraday rotation due to interactions with background magnetic fields, and (c) interactions with 
pseudoscalar fields [22|. The interactions with foregrounds leads to a frequency-dependent effect. The same is true of 
Faraday rotation, where the frequency dependence is ~ [H, [13, H^, [1^ , while the interactions with pseudo-scalar 
fields are frequency-independent. These rotations change the CMB Stokes parameters as 

(Q±iC/)(n) e-*"(")(g±ii7)(n) => S{Q ± iU){n) = ~ia{h){Q ± ilJ){n) , (21) 

where Q;(n) is the integrated rotation of the plane of polarization along the line-of-sight n. Comparing Eqn. (j2ip with 
Eqn. U) it is clear that the cosmological rotation is captured in our rotation distortion field ^(n). Currently there 
are no direct constraints on spatially varying rotation, however indirect constraints from the upper limit of B-modes 
exist [13] ■ The constant rotation has been constrained to be less than few degrees [l3, [13 • From the QUAD data [1H| 
the constant rotation has been constrained to 0.55 ± 0.82 (random)±0.5 (systematic) degree. Planck can reach the 
sensivity of O.ldeg, while futuristic experiment like CMBPol will be able to detect rotation angle as small as 0.005 
degree 

Patchy reionization: Inhomogeneous reionization introduces several effects in the CMB map. First, ionized materials 
scatter the local CMB temperature quadrupole, thus generating extra polarization on large-scales [s^, [Slj. Second, 
extra temperature and polarization anisotropics are generated from peculiar motion of ionized regions during the 
entire reionization process, introducing a kinetic Sunyaev-Zeldovich signal [s^, [ssl . [s^ . [ssl [sgI . Wf\ . Finally, the 
primordial anisotropics are suppressed by a factor of e~'^'^"\ where T(n) is the line-of-sight dependent optical depth 
to recombination. Scattering which mixes photons from different directions, is not a line-of-sight distortion and is not 
modeled in Eqn. ([1]), while the screening effect from patchy reionization is a line-of-sight effect and is indeed captured. 
The observed polarization Stokes parameters due to the screening effect are 

(Q±iC/)(n) = e-^^^\Q±iU){ii)=>5{Q±iU){n)^-T{h){Q±iU){h). (22) 

Comparing the above equation with Eqn. ([T]) shows that the screening effect from patchy reionization is captured by 
our modulation parameter a(n). For typical patchy reionization models, i5T(fi) is of the order of 10"'^ and hence it 
will be challenging to reconstruct it. The corresponding B-modes arc about an order of magnitude smaller than the 
Icnsing induced B-modes [l^, [11] . 
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Lensing: The path of the CMB photons get deflected by the intervening matter between us and the last-scattering 
surface. The remapping of CMB photons due to this so caUed weak gravitational lensing of the CMB can be written 

as filllilii 



{Q±iU){h) (Q±iC7)(n + V$ + V X fi), (23) 
which after Taylor expansion can be written as 

6{Q ± iU){n) = V$ • V(Q ± iU){n) + V xQ- V(Q ± iU){h) + ... . (24) 
Here $ is a weighted projection of the gravitational potential along the line-of-sight and is written as 

$ = -2 1 dD ^^*^~^^ ^{Dh, D) , (25) 

where is the comoving distance to the last-scattering-surface. The gradient part of lensing, V$, is generated by 
the linear order density perturbations, while the curl part, V x 57, comes from the second order density perturbations. 
The curl part is much smaller than the gradient part and is usually ignored. The form of Eqn. (j24[) when compared 
with Eqn. ([1]) shows that the gradient type of lensing is captured in our deflection field p2 (n) , while the curl part of 
lensing is same as our deflection distortion pi (fi) . 

Other cosmological distortions: Although there are several other CMB distortions in the literature such as spatial 
variation of the fine structure constant [4l| and vorticity from topological defects [13], we will not discuss their cor- 
respondence here as our objective here is to just point out that there are several cosmological distortions that are 
captured in Eqn. (1). We do want to emphasize that not all the distortions in Eqn. (1) correspond to cosmological 
distortions. For example as far as we know there are no known physical processes which can covert CMB temper- 
ature to polarization along the line-of-sight, hence there is no physical analog of temperature leakage distortions 
(711 72), (c^i, 0^2), and q. However we will now see that instrumental systematics can generate all the 11 distortion 
types, including leakage from temperature to polarization. 



B. Instrumental Systematics 

Time varying errors produced by detector systems translate into position dependent distortions in the map. Even 
if the detector response is not time varying, for multiple detector systems if response of the different detectors have 
different systematics then that also results in position dependent distortion in the map. A map making process can 
be encapsulated as 

d = As + n, (26) 

here d is a vector containing time ordered data; A is a pointing matrix and encodes sky pixel location for each 
observation, scanning strategy, and beam effects; s is a map we want to reconstruct and fi is instrument noise. The 
relationship between the instrumental systematics and map distortion depends on scanning strategy. For a given 
experiment design and scanning strategy, starting from detector error models, one can calculate the distortions in the 
map. 

Now we point out correspondence between 11 distortion fields in Eqn. (1) and various instrumental effects (for 
details refer to [l^)- The effect of detector gain response miscalibration is captured by our modulation field a(n). 
The effect of errors in the detector orientation is described by the u!{n) field. If the two polarization states are 
measured separately, then the differential gain is captured in our monopole leakage field (71 (n), 72(11)). Apart from 
the errors in the detectors there are also errors due to finite size of the beam. Let us consider an example of beam of 
the form 

^' = W(T^ i - 2^ [iTTef- + -JT^J \ ' 

where b is the beam offset, a is the mean beam-width, and e is the ellipticity. These parameters could be different for 
the different polarizations, and only the difference in the beams enters into the CMB measurements. Assuming any 
distortion of the beam are relatively small compare to typical scales of the beam on either direction, one can show 
the correspondence of beam distortions with the distortion parameter's considered in the paper 

(TP = (b„ +b6)/2, 
ad = (b„-bb)/2, 
es = (ea + eb)/2 , 

g = (ea - eb)/2 , (28) 
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FIG. 1: Estimator Gaussian noise N'^g{L) as a function of multipole L for the EB estimator. Upper panels show the noise for 
modulation a, rotation co, monopole leakage (71, 72), and spin-flip (/i, f^); while the lower panels show the noise for deflection 
p, dipole leakage (di,d2), and quadrupole leakage q. For the left panels we assume the instrument with polarization noise 
Ap = 9.6/iK-arcmin and beam FWHM = (81n2)^/^a = 8'. For the right panel we assume a CMBPol-like instrument with 
polarization noise Ap = f .4^K-arcmin and beam FWHM = (81n2)i/V = 4'. Calibration a is dimensionless, rotation lu is in 
radians, deflection p and dipole leakage (di, d2) are in the units of Gaussian beam a; and quadrupole leakage q is in the units 
of . For reference we also show the Gaussian noise for lensing deflection estimator. 



Hence all the 11 distortion fields in Eqn. (1) correspond to some sort of instrumental systematics [1. 



V. RESULTS AND DISCUSSION 

Here we discuss the prospects of using our estimators for detecting the distortions in the CMB. Wc study the 
estimator variance for two experimental setups, (1) experiment with noise sensitivity of Ap — 9.6^K-arcmin and 
beam FWHM a ~ 8', typical of upcoming ground and balloon based CMB polarization experiments (EXP-balloon 
hereafter) and (2) a futuristic CMBPol like satellite experiment with noise sensitivity Ap = 1.4/iK-arcmin and beam 
FWHM (7 = 4' (CMBPol hereafter). Wc choose a standard fiducial model with a flat AC DAI cosmology, with 
parameters described by the best fit to WMAP5 [l3|, given by Qb ^ 0.045, Qc = 0.23, Ho = 70.5, = 0.96, = 0.0, 
and T = 0.08. We calculate the theoretical lensed and unlensed CMB power spectrum from publicly available code 
CAMB'^ [4^. We assume zero distortions as our fiducial model, except for the lensing of CMB by the large-scale 
structure whose effect we include in the variance by using the lensed CMB power spectrum. For the experiments in 



^ http://camb.info/ 
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FIG. 2: Total signal-to-noise S/N for the EB estimator as a function of maximum multipole L,nax- We assume that the 
distortion field power spectrum is white noise above a coherence scale aj}, with rms amplitude A-d- The exact form of the 
distortion power spectrum is given by Eqn. (|19|l . The left panel assumes the instrumental noise and beam properties of EXP- 
balloon experiment, while the right panel assumes that of CMBPol experiment. In both the panels we have assumes full-sky 
coverage. 
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0.045 


Dipole leakage di 
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5.23 
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0.0077 


0.0053 


Dipole leakage d2 


1.7 


5.38 




3.8 


11.11 


0.0077 


0.0056 


Quadrupole q 


1.8 


4.1 




3.32 


3.55 


0.0124 


0.0394 


Deflection pi 


38.2 


19.4 




132.7 


40.3 


0.75 


0.53 


Deflection p2 


4.4 


15.5 




10.8 


24.8 


0.098 


0.57 



TABLE II: Here we show that using the EB correlations, the distortions will be detected at high significance before the distortions 
show up as B-modes. In column 1 we show the distortion type in consideration. For a given distortion type, the numbers in 
column 3 show the maximum allowed distortion rms A^""^ above which the distortions will be detected as B-modes with an 
experiment which is sensitive to r = 0.005. These numbers are from [l3|. In column 2 we compare the minimum detectable 
distortion rms amplitude A^'" (for experiments EXP-balloon and CMBPol) using our EB estimator with the maximum allowed 
distortion rms ^4^^" presented in column 3. The minimum detectable rms Ag'" was obtained using Eqn. (|20p and seeking 
the value of At> for which total signal-to-noise becomes unity. We have considered two choices of coherence length a-p = 10' 
and a-D = 120'. For same a-p, the numbers in column 2 are much larger than the numbers in column 3, this implies that the 
distortions will be detected using EB estimator before they will be detected as B-mode power spectrum. 

consideration, the EB estimator is most sensitive so we will restrict our discussion to that. 

In Fig. Il]) we show the Gaussian noise N^^{L) for all the 11 distortion fields. Notice that noise N^^{L) for 
distortions a, q, and {pi,p2) whose filter /ebO-i, I2) contains sin(iy9ij — (ySi^) scales as N{L) oc for low L. This scaling 
can be understood by noting that N{L) oc ^.^^j"'^ and for small L, li « I2 and hence sin((^ij, (^i^) « V'lib ^ L/li. In 
a similar way the sealing of other distortions noise can be understood, for example noise N{L) for rotation u whose 
filter /£;s(li,l2) contains cos{ipi-^ , (pi^) is roughly constant N{L) w constant at low L. This is because for small L, 
h w I2 and hence cos((/'ij, (/Sij) w 1. 

For a given multipole i, if the distortion power spectrum is larger than the Gaussian noise N'^ {£) then the estimator 
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can reconstruct the distortion map with (S/N)i > 1. However if the distortion power spectrum is smaller than the 
Gaussian noise than the distortion map cannot be reconstructed but one can still statistically detect and constrain its 
size. Later in the section we consider an application of the estimator for doing diagnostics for the primordial B-modes. 
There we will quantify the level of distortions that can be detected using our estimators, and how they compare to 
other methods. 

In Fig. [2] we show total signal-to- noise for all the 11 distortions as a function of £max for EXP-balloon and CMBPol 
experiments. To calculate the S/N we have assumed that the power spectrum of distortion is given by Eqn. (jl9p . 
The S/N scales as A^. In the figure we have chosen rms At> = VO.Ol for EXP-balloon experiment and Ad = VO.OOl 
for CMBPol experiment^. The coherence scale is taken to be cr-p ~ 10' for both the experiments. For same distortion 
rms, among all the distortions, monopolc leakage has the highest S/N. This is because the amplitude of temperature 
fluctuations is much larger than the polarization amplitudes and monopole leakage represent a distortion where 
temperatures gets converted into polarization. 

Although the estimators are unbiased when averaged over CMB realizations, for a given CMB realization (which is 
the case in reality) the reconstruction is noisy even for the cosmic variance limited experiment with perfect beam. In 
Fig. ([3]) we show the reconstruction of distortion for a given CMB realization. We consider a modulation distortion 
and show the reconstruction ability for EXP-balloon like experiment. The simulation patch size is 30 x 30 square 
degrees. We assume the modulation to be a Gaussian field with the underlying power spectrum given by Eqn. (fTOj) . 
We have chosen the coherence cr-p and rms A-p in such a way that the modulation induced B-modes resemble the level 
and shape of primordial B-modes that EXP-balloon is sensitive too. In the lower left panel we compare the B-modes 
generated from the modulation field with the primordial B-modes with r = 0.01 and r = 0.1; and lensing induced 
B-modes. The upper left map shows the realization of modulation field a'"''"*(fi) that we use for the simulation. 
The upper right panel shows the reconstructed modulation field for a single CMB realization. It is clear that the 
reconstruction is noisy however the reconstructed modulation traces the features of input modulation. We have 
smoothed the input and reconstructed maps with Gaussian beam FWHM of 1° so that it is easy to see the common 
features. The middle panel shows the CMB E and B fields which were used to reconstruct the modulation field. The 
middle left map shows the E-polarization of the CMB after being distorted using input modulation field a'"^"*(n). 
The middle right map shows the B-mode polarization generated by the modulation field a'"^"*(n), note that these 
B-modes are solely due modulation of primordial E field as we have assumed zero primordial B-modes for our fiducial 
model. In the lower right figure we compare the reconstructed binned power spectrum with the power spectrum of 
input modulation a'"P"*(n). 



Primordial B-mode Diagnostics 

As we discussed in the introduction, detection of primordial gravitational wave provides information about the 
early universe. So far primordial gravitational waves have not been detected, there are only upper limits on their 
amplitude, see (4^ . for a current observational bounds on the level for primordial gravitational waves. If the B-modes 
are observed what diagnostics test can be done to make sure that they are primordial [i^? As we have already shown 
that many distortions in the CMB which produces B-modes also produce off-diagonal correlations in the CMB. Using 
the off-diagonal correlations in the CMB, we have also constructed estimator to reconstruct the distortions. Now wc 
show that our estimators are more sensitive to distortion than the B-modes and any non-primordial B-modcs which 
can be detected, can be detected using our estimators at much higher significance. 

In Table [III we compare the minimum detectable level of distortion rms using our EB estimator with the 

maximum allowed distortion rms for primordial B-modes detection. We consider two experiments EXP-balloon and 
CMBPol, which are sensitive to scalar-to-tensor ratio of r = 0.05 and r = 0.01, respectively. Wc find that the EB 
estimator will be able to see any spurious contamination before it gets detected as B-modes. In column 1 we tell the 
distortion type in consideration. For each distortion type, the column 3 shows the maximum distortion rms Ap^'^ 
above which the distortion will be detected as B-modes for an experiment that is sensitive to r = 0.005. In column 
2 we show the ratio of maximum allowed rms given in column 3 and the minimum detectable rms using the our EB 
estimator. For EXP-balloon, which is typical of upcoming balloon based experiment, we assume sky fraction to be 
1%. For CMBPol experiment we assume full sky fraction to be fsky = 0.8. For distortion field we have considered 
two choices of coherence length a-p- For a given coherence ax>, the fact that numbers in column 2 are larger than 



* The values of rms in Fig. [2] wore mostly chosen for visual purpose. Those rms values roughly correspond to 10 a detection of deflection 
distortion p2. 



10 



Input 



Reconstructed 




0.2 -0. 




I'm m ■ ■* 




"T 1 r 



10 X 




m 

< 



0.1 



0.01 



0.001 



Primordial (r=0.01) 
Primordial (r=0.1) 
A_=9.6 |iK-arcmin Modulation 
a=8' Lansing... 
o,=120', A,=0.1 




0.1 



0.01 



O 



0.001 r 



0.0001 



10 



100 



1000 



1e-05 



fsky=0.02 
a=8' 

Ap=9.6 )iK-arcmin 



Input 

Reconstructed :----X- 




10 



100 



FIG. 3: The upper left map shows the input modulation map a'"''"*(n), the upper right map shows the reconstructed a(n) map, 
the middle left map shows the E map after primordial E map has been distorted via input modulation map (upper left map) , the 
middle right map shows the B-polarization field generated by input modulation. In our simulations we have assumed primordial 
B-modes to be zero. The lower left map compares the B-modes generated by modulation with the primordial B-modes (for 
r — 0.1, r — 0.01) and lensing induced B-modes. The lower right panel compares the power spectrum of input modulation field 
with the reconstructed binned power spectrum. The size of the simulation patch was 30 x 30 square degree. The instrumental 
noise Ap = 9.6/iK-arcmin, and beam FWHM cr = 8' correspond to EXP-balloon, an experiment representative of upcoming 
balloon and ground based experiments. 
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FIG. 4: Ratio of tensor-to-scalar uncertainty Ar without and with cleaning as a function of distortion rms A-p- We have 
assumed convergence length to be au = 10'. The left panel is for the EXP-balloon experiment, and right panel is for the 
CMBPol experiment. The ratio of uncertainty is independent of tensor perturbations, and sky fraction. 



unity says that the distortions will be detected at higher significance using EB correlations than using B-niode power 
spectrum of the CMB. 

Detection of distortion V through these distortion estimators clearly determines the presence of iron-primordial B- 
modes. In addition it may mean two things 1) either the source of the distortion is instrumental systematics or 2) the 
distortion signal is cosmological (although not primordial^). If the source of the detection is instrumental systematics, 
the significance of the detection will drop as the instrumental systematics is controlled better (or another instrument 
with better systematics (corresponding to distortion V) control will not see the detection using the V distortion 
estimator^). However if the systematics is controlled well below the noise level N'^{L) of V distortion estimator, and 
the signal is still present, then the detection should be attributed to cosmological signal. All cosmological distortions 
in the CMB probe interesting cosmological processes and in some cases probe physics beyond the standard model [3 ■ 
Hence detection of cosmological distortion will be extremely exciting. 

So far we have only shown that the distortions can be detected in the power spectrum. However if the S/N per 
mode is much larger than unity, our estimators can be used to reconstruct the maps of the distortion fields. The 
reconstructed distortion then can be used to clean the CMB B-modes by subtracting the spurious B-modes from the 
observed B-modes, hence improving the sensitivity for the tensor B-modes. Fig. |4] demonstrates the improvement in 
the uncertainty of tensor-to-scalar ratio Ar due to the cleaning. We show the ratio of tensor-to-scalar uncertainty Ar 
without and with cleaning as a function of the distortion rms A-p. In Appendix B we describe the B-modes cleaning 
process. While both Ar^ith cleaning and Arno cleaning depend on the value of tensor B-mode power spectrum and sky 
fraction fsky, the ratio is independent of them. It is clear from the figure that larger the rms of the distortions, larger 
the reduction in the error of tensor-to-scalar Ar due to cleaning ^. 



VI. SUMMARY 



We have considered general line-of-sight distortions in the CMB. These distortions can be described by 11 fields 
and their effect on CMB polarization is given by Eqn. (1). We show that these fields correspond to various known 
instrumental systematics and cosmological signals. Example of cosmological signals include rotation of plane of 
polarization due to magnetic fields or parity violating physics, screening effects of the patchy reionization, lensing 
of CMB by the large-scale structure. Instrumental systematics include gain fluctuation, differential gain, pointing, 
beam ellipticity, differential beam ellipticity, detector rotation. All these distortions generate B-modes of CMB and 



^ Primordial fluctuations do not generate (EB) correlations, and hence our quadratic estimator will not see anything. 

® The Distortions due to systematics will not correlate between different experiments. 
Note that if the rms is below the requirement from tensor B-mode detection (i.e. column 3 of Table II), one docs not reduce the 
uncertainty in r by much from the cleaning process, for example for CMBPol. For example for CMBPol the the uncertainty Ar reduces 
by only factor of 5 for the deflection pi and is even smaller for other distortions. 
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also generate unique ofF-diagonal correlations in the E and B fields. The way the different modes of E and B get 
correlated depends on the distorting field V and can be compactly written as (£^(li)-B(l2)) = /-^(li, l2)^(li + 12)- 
We construct quadratic estimators for all the 11 distortion fields which utilizes the unique correlation signature in 
the observed CMB fields to reconstruct the distortion fields. We compute the variance of the estimators and show 
the level of distortions that can be detected in upcoming and future experiments such as CMBPol. As an example of 
our estimator we show that our estimator can detect the distortions before they can be detected as spurious B-modes 
in the power spectrum. This is particularly important diagnostics for primordial B-modes, since distortions induced 
B-modes act as contamination for primordial B-modes detection. Since there are only 11 distortions, it is worthwhile 
to look for them in CMB as systematics check or if instrumental systcmatics is well controlled, to look for cosmological 
signals. 
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APPENDIX A: ESTIMATING MULTIPLE DISTORTIONS SIMULTANEOUSLY 

111 the main text we discussed estimator for different distortion fields, while estimating them one at a time i.e. in 
the absence of any other distorting field. In this appendix we discuss how can we estimate multiple distortion fields 
simultaneously. For simplicity we will restrict the discussion to the EB estimator, the generalization to include other 
CMB fields is straight forward. To generalize our EB estimator to multiple field we start by defining a N-hy-N 
matrix FJ^'^ for each multipole £, 

= l^fUh,£2){C-X''.f^B{h,i2)(C-X'', (Al) 

where V and V run over all the 11 distortions and lensing, {a, lu, 71, 72, /i, /2, c?i, rf2, Q,Pi,P2, 0}- An estimator in the 
presence of single distorting field was unbiased, 

(P(L))^,,^ = V{L) . (A2) 

In the presence of multiple distorting fields the estimator as constructed in Eqn. (10), may not be unbiased. In general 
the Eqn. (|A2p gets modified as 

{'D{L))cMB = V{L) + ^ ^ , (A3) 

where T) is an estimate of T) in the presence of multipole distortions effect in the CMB. The level of bias depends on 
the off-diagonal Fisher matrix Fj^^ and the amplitude of distortion fields V . Eqn. (jA3P can be inverted iteratively 
to get an unbiased estimate of a field D. The number of iterations required to converge depends on the signal levels 
and correlation between the fields. The correlation between two fields at a give multipole L can be calculated as 

Cv-D' = , ' (A4) 

In Table UlTl we give the Fisher matrix at L = 100 for EXP-balloon experiment. Most of the distortion fields are 
uncorrelated^, however there are some fields which arc highly correlated. For example the largest correlation is 
between quadrupole leakage q and dipole leakage c?2j with the correlation coefficient Cqd^ — .9. The second largest 
correlation is between deflection pi and rotation lo, with the correlation coefficient Cp^^ = .8. The next three largest 
are Cp^a = —0.2, C^^j = .15 and Ca-fi ~ .14. 

The variance of the distortion T) will be given by (^F^^^^^. Multiple field formalism for other estimators 
(TT, TE, EE, TB) and the minimum variance estimator is a straight forward generalization of the EB estimator 
we discuss here. 



As we have shown in Sec. IIVI lensing deflection is same as deflection p2. Hence as expected the correlation coeflicicnt between lensing 
and p2 is unity, Cp^Xens.d = 1- 
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TABLE III: Fisher matrix Ff ^' at £ ^ 100 for EXP-balloon experiment with Ap = 9.6/iK-arcmin and beam FWHM= 8'. 



APPENDIX B: SELF CALIBRATING THE CMB 



In the main text we have only shown that the non-primordial contribution to the B-modes can be detected using 
the EB estimators. Here we turn to the question that whether we can clean the CMB from the distortions and 
hence improve the sensitivity to primordial B-modes. It is also useful to correct the systematics induced distortions 
as they also affect lensing reconstruction from CMB [2l|, |4^ . If the estimator S/N per mode is greater than unity, 
then the estimators can be used to reconstruct the map of the distortion fields without any assumptions about the 
statistics of the distortion fields (i.e. The distortion field does not have to be Gaussian). Once the distortion fields 
are reconstructed one can clean the observed CMB by un-doing the distortions in CMB. The process is similar to 
de-lensing [3§, [13, IH, SI . 

In principle one can clean the CMB field from distortions iteratively. For each iteration, the first step is estimating 
the distortion filed V using the observed CMB E and B fields, i.e. using Eqn. [TOl Then constructing an estimate of 
distortion induced B-modcs given the distortion field 2?, and observed CMB field X can be calculated as 



B{1) 



(27r)2 



W(ll,l2) 



(Bl) 



where X = E,Y = B ior V = a,uj, (/i,/2), {pi,P2) and X = T, F = S for P = (71,72), (^1,^2), 9- The windows 
are given in Table [J The power spectrum Cf^ {cleaned) of cleaned B-modes Bdeanedi^) = ^(1) ~ ^(1) is thus given 
by 
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We can repeat the above process of cleaning by using the cleaned Cf^ {cleaned) given by Eqn. (jB2[) as the new 
observed Cf^ and repeating the cleaning process until Cf^ {cleaned) has converged. In each iteration of cleaning 
process the Gaussian noise of the estimator Nf^ reduces until it converges to maximum likelihood estimator noise. 
The decrease in the estimator Gaussian noise can be understood by noting that the B-modc power spectrum acts as 
noise source for the quadratic estimator. By iteratively cleaning the B-modes, we are reducing the noise source for 
the estimator until B-modes are nothing but primordial (i.e. tensor) B-modes. 
The uncertainty is the tensor-to-scalar Ar is given by 
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where Arno cleaning is the uncertainty in r when no cleaning was performed, while Ar^ith cleaning is the uncertainty 
after the distortions have been cleaned. Note that the ratio of uncertainties with and without cleaning is independent 
of tensor perturbations, and sky fraction covered. 



